
% 基于TimeSpecialist man探究环境因素对时间的影响
clc,clear,close all
load("map_tokyo.mat")
%% 已知参数
Humidity = 0.59;
hight = y_altitude;
theta = 0; %风与运动员垂直方向的夹角,单位是度
V_wind = 3.*ones(1,length(GR));
Temperature = 288.75;
tensity = 500;
CP = 411.1;
W = 20720;
% figure()
% plot(GR)

%%  V_wind变化, theta=0
s = 1:1:4;
for n=1:length(s)
V_wind = s(n).*ones(1,length(GR));
[P_plan,P_real,V,delta_T,fatigue,t(1,n)] = get_P(GR,CP,W,tensity,hight,V_wind,theta,Humidity,Temperature,turn_sign);
%需要画图可以去掉相应注释
% title('GR-x')
% figure()
% plot(V)
% title('V-x')
% figure()
% plot(P_real)
% title('P-x')
end
figure()
plot(s,t,'o-')
xlabel('V_wind')
ylabel('t')
title('TimeSpecialist的时间t随V wind的变化')
%% theta变化, V_wind=3
s = 0:2.5:7.5;
for n=1:length(s)
theta_Vary = s(n);
[P_plan,P_real,V,delta_T,fatigue,t(1,n)] = get_P(GR,CP,W,tensity,hight,V_wind,theta_Vary,Humidity,Temperature,turn_sign);
end
figure()
plot(s,t,'o-')
xlabel('theta')
ylabel('t')
title('TimeSpecialist的时间t随theta的变化')
%% Humidity变化
s = 0.57:0.01:0.60;
for n=1:length(s)
Humidity_Vary = s(n);
[P_plan,P_real,V,delta_T,fatigue,t(1,n)] = get_P(GR,CP,W,tensity,hight,V_wind,theta,Humidity_Vary,Temperature,turn_sign);
end
figure()
plot(s,t,'o-')
xlabel('Humidity')
ylabel('t')
title('TimeSpecialist的时间t随Humidity的变化')
%% Temperature变化
s = 286.75:1:289.75;
for n=1:length(s)
Temperature_Vary = s(n);
[P_plan,P_real,V,delta_T,fatigue,t(1,n)] = get_P(GR,CP,W,tensity,hight,V_wind,theta,Humidity,Temperature_Vary,turn_sign);
end
figure()
plot(s,t,'o-')
xlabel('Temperature')
ylabel('t')
title('TimeSpecialist的时间t随Temperature的变化')
